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Abstract 



We investigate cosmological evolutions of the bulk scalar field <p(t) and the radion d(t) 
in five-dimensional dilatonic two branes model. The bulk potential for the scalar field 
is taken as the exponential function Vbulk °c exp(— 2\/26<^), where b is the parameter 
of the theory. This model includes Randall-Sundrum model (with 6 = 0) and five- 
dimensional Hofava-Witten theory (with 6 = 1). We consider matter on both branes 
and arbitrary potentials on the branes and in the bulk. These matter and potentials 
induce the cosmological expansion of the brane as well as the time evolution of the 
bulk scalar field and the radion. Starting with full five-dimensional equations, we derive 
four-dimensional effective equations which govern the low-energy dynamics of brane 
worlds. A correspondent five-dimensional geometry is also obtained. The effective four- 
dimensional theory on a positive tension brane is described by bi-scalar tensor theory. 
If the radion is stabilized, the effective theory becomes Brans-Dicke (BD) theory with 
BD parameter 1/26 2 . On the other hand, if the scalar field is stabilized, the effective 
theory becomes scalar-tensor theory with BD parameter 2 (3$+i) i-y(t) wnere V 9 is the 
BD field defined by radion d(t). If we do not introduce the stabilization mechanism for 
these moduli fields, the acceptable late time cosmology can be realized only if the dilaton 
coupling 6 is small (6 2 < 1.6 x 10~ 4 ) and the negative tension brane is sufficiently away 
from the positive tension brane. We also construct several models for inflationary brane 
worlds driven by potentials on the brane and in the bulk. 
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1 Introduction 



String Theory and/or M-Theory are thought to be the most promising candidate for "The Theory 
of Everything". One of the interesting and embarrassing features of such theories is that they 
are consistently formulated in more than four-dimensions. On the other hand, at least at low 
energy scales, our spacetime seems to be four-dimensional. Hence we need some compactification 
mechanisms of extra dimensions in order to retrieve our four-dimensional spacetime at low energy 
scales. 

Of all compactification mechanisms of the extra dimensions, Hofava-Witten (HW) theory is a 
fascinating onepQ[2]. Hofava and Witten showed that the eleven-dimensional limit of M-theory 
compactified on S 1 /Z2 orbifold can be identified as the strongly coupled ten-dimensional E% x Eg 
heterotic string theory with two orbifold fixed planes. Witten also showed that there exists a 
consistent compactification of this theory where the universe appears as a five-dimensional spacetime 
with two boundary branes ,3|. Then the five dimensional effective action for this theory was derived 
[I] and the cosmological evolution of the branes has been studied |S]-jHl- 3 In this model, the bulk 
scalar field is necessarily associated with the compactification. 

The bulk scalar field also plays an important role in phenomenological models for brane worlds. 
For example, the bulk inflaton models [H]-J2]> the self-tuning models for cosmological constant 
^H]- JZj and the dilatonic thick domain wall models ^B] -|2J| have been studied as the extension of 
Randall- Sundrum (RS) models |23| . Cosmological models with bulk scalars have been actively 
considered |24| . |25| . In the context of RS model, stabilization mechanism of the distance between 
two branes using bulk scalars has been investigated |2l?]-|33|. 

As well as the bulk scalar field, there is another important scalar degree of freedom in two branes 

3 Note that the time evolution of their cosmological domain wall solution is driven by the bulk scalar field, not by the 
matter on the brane. So their solution is a "vacuum cosmological solution", which is different from the cosmological 
solution in our usual sense. 
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model, that is, a radion which describes the inter-brane distance. The radion plays a crucial role in 
gravitational interactions on the brane. For example, in RS model, the radion acts as a Brans-Dicke 
(BD) field and the observational constraints on BD parameter gives strong constraint on the model. 
In RS model, the dynamics of the radion has been eagerly studied |31]-|4L)|. 

The purpose of this paper is to investigate the cosmological consequences of the bulk scalar field 
and the radion. We consider a model where the bulk potential for the bulk scalar field <j) is given 
by an exponential function V^ik oc exp(— 2y/2b(p). This model includes HW theory (with 6 = 1) 
and RS model (with 6 = 0). We will consider matter on both branes and arbitrary potentials on 
the branes and in the bulk. These matter and potentials induce the cosmological expansion of the 
brane as well as the time evolution of the bulk scalar field and the radion. It is well known that 
in many higher-dimensional theories, these "moduli fields" are potentially dangerous because time 
evolution of moduli fields leads to the time variation of fundamental constants. Then it is important 
to know how serious is the moduli problem in this model. We also consider the phenomenological 
implication of the potentials mainly focusing on the inflationary scenario. 

In order to know the cosmological consequences of the bulk scalar field and radion, we should 
derive the effective equations which govern the cosmological dynamics of these fields. One approach 
is to construct an 4D Effective action by the integration of the action with respect to the extra 
dimension 0^-01]. A "moduli space" approximation is an example for such an approach. In this 
approximation, the constant moduli parameters of the static solution is replaced by time-dependent 
functions and integration over the extra dimensions gives four-dimensional effective action. But, 
it is known that the moduli space approximation may not always be consistent |45j-|47j and the 
integration over extra-dimension must be treated with special caution 48 . 

Hence, in this paper, we will derive four-dimensional effective equations for the radion and the 
bulk scalar field starting with full five-dimensional equations. In order to do so, we use the low- 
energy expansion scheme. The low-energy means that the energy density of the matter on the 
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brane is much lower than the tension of the brane. This relation holds at late times in cosmological 
evolution. Under such an approximation, we can derive the effective four-dimensional equations for 
scale factor, radion and bulk scalar field. We can also determine the correspondent five-dimensional 
geometry. 

The organization of this paper as follows. In Sec. II, we review the setup of our model and 
introduce the low-energy approximation. Using this approximation, we give the background solu- 
tions. In Sec. Ill, we solve the first-order equations and derive the solutions for full five-dimensional 
metric at this order. Using these solutions, we show that the equations for scale factor on the brane, 
bulk scalar and radion can be derived. Sec. IV is the main part of this paper. Effective equations 
which govern the low energy dynamics of the two branes system are presented. The comparison 
with scalar-tensor theory is performed and the observational constraints are obtained. Sec. V, 
several inflationary scenarios are discussed. Sec. VI is devoted to conclusions. 

2 Setup 

2.1 Action and Equations of Motion 

We consider the bulk action of the form 



Jm 5 

where k 2 denotes the five-dimensional gravitational constant and a denotes the tension of the 3- 
brane. The bulk potential is taken as 



Here we split Vbulk i n t° an explicit, exponential potential for <f> and an arbitrary potential V{4>) 
according to Ref. |49j . The parameter b is determined by the theory. If we choose 6 = and 
V((f>) = 0, we retrieve the action for the Randall-Sundrum model and if b = 1 and V{(p) = 0, we 
retrieve the action for Hofava-Witten model, which were derived in [I]. 
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We consider a five-dimensional spacetime with structure M5 = S 1 /Z2 x M4 where S 1 /Z% is an 
orbifold and M4 is a smooth manifold. We use the coordinate y to parameterize the orbifold and 
the range of y is chosen as — r < y < r. Because this spacetime has Z2-symmetry (i.e. we identify y 
with — y), there are two endpoints in this spacetime. These two endpoints are called orbifold fixed 
planes, which are 3-branes. A positive tension brane is set at y = and a negative tension brane is 
set at y = r. We take the actions for them as 

Sbrane = S+ + S- 

S+ = / ^d 4 Xy/^g^U brane + / d A x^f^giC m , (3) 

JM^> JM^> 

S- = d 4 Xy/=g^U brane + / d 4 Xyf^giC m , 

JMf> JMl 2 ' 

where the potentials are given by 

Ubrane = ~e~^<T - U(<f>), 



V2 



U brane = ^- e -^a-U(4>). (4) 

The action for the positive tension brane and the negative tension brane is denoted by S+ and 
S-, respectively. U{<f>) and U{4>) are arbitrary potentials. C m and C m are the Lagrangian density 
of the matter on each brane. We assume that they do not couple to the bulk scalar, that is, 
5 C m / 5(f> = 5C m /5(p = 0. It is true that in HW model, there are field strengths of Eg and Eq gauge 
field which couple to the bulk scalar field, but in this paper, we replace them with the energy density 
of perfect fluids . 

Using the action, we derive the Einstein equation 

n M Jl(T M , rp M \ /K\ 

u N — K [J-bulk N ' 1 brane NJ- \°) 



It is also interesting to consider the coupling of the gauge field and the scalar field. It is left for the future work. 
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The energy-momentum tensor can be divided into two parts; 



r M 
- bulk N 



rp M 

brane N 



1 




v- 


-94 , 


v- 


~9b 


+ 


V2 

^2 



V M </>V M </> + ( fo 2 - \ ) e^V - F(<£) 



V2 



^ dia0(O, a, a, a, a) + c&a^O, -p - U,p - U,p - U,p - U) 



diag(0, a, a, a, a) + diag(0, —p — U,p — U,p — U,p — U) 



S(y-r) YB) 



The energy density and the pressure of the perfect fluid on the positive tension brane are denoted 
by p and p and those on the negative tension brane are denoted by p and p. The equation of motion 
for the scalar field is given by 



vv + v&(*»-§) e -*»v + i*: 



-2be~^a + k 2 ^J 8{y) + 2be~^a + k 2< ^\ 6(y - r) 



V^95 

Now we take the five-dimensional metric as 

ds 2 = e 2 ^dy 2 - e 2 ^dt 2 + e 2a ^ '5 l3 dx l dxK 



(7) 



Einstein equations and the scalar field equation in terms of a, (3, 7 and <p are presented in Appendix 
A-l. We expand the functions around y = and y = r as 



a(y,t) = a W(t)+a^(t)\y\ + ^ 2 Ht)y 2 + ..., 
a(y,t) = a^(t)+u^(t)\r-y\ + \aW(t)(r-y? + 



(8) 
(9) 



From above expansions, we can read that a" gives 25{y) — 28(y — r). Then from Einstein equations 
and dS^ can be written in terms of p, p, U and p, p, U respectively. Then we can get the 
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junction conditions; 



a«(t) 
5«(t) 



l r 

6 

1 

6 L 



+ K (p+U) 



V2ae-^ - K 2 {3p + 2p-U) 



y=0 



bae -V2b<f> _ ± K 2 

2 



i „dcr 



y=0 



1 r 

6 

l r 

6 



V2ae~^-K 2 (p + U) 
yfae-y/M* + K 2 (3p + 2,5 - £/) 



1 dU 
2 



(10) 

(11) 
(12) 

(13) 

(14) 

(15) 



y= r 



2.2 Low Energy Expansion Scheme 

In this subsection, we formulate the low energy expansion scheme to solve the equations 5 . At late 
times, the energy density and the pressure of the perfect fluid is much smaller than the tension of 
the brane. Hence we can impose the following conditions, 



a > K 2 p, K 2 p, cr~ l V, k 2 U, n 2 U. 



Then, we can define a small parameter of the system as follows, 



K 2 p 
a 



(16) 



(17) 



Thus for £ < 1, we can solve five-dimensional Einstein equations perturbatively by expanding a 
metric and scalar field as 



a(y,t) = ac(p)(y,t) +a (1) (y,t) +a (2 )(y,i) H 



(18) 



where atr^ is the function of the order 0(e l ). Note that the lower indices au\ are different from the 

upper indices a'*' which have been defined in eqs.(jSJ), ©• We expand f3(y,t), 7(2/, t), and 4>(y,t) 
Generalizations to covariant form of this iteration scheme are performed in |50| . |51| 
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in a similar way. The important point is that i-derivative of the lower order solution is the same 
order as y-derivative of the higher order solution; 

<%-i) ~ (19) 

This is because the derivative with respect to y is of the order a' ~ a for e <C 1 from the junction 
conditions. On the other hand, from the Friedmann equation, we expect that the derivative with 
respect to t is of the order a 2 ~ K 2 ap. Then the derivative with respect to t is suppressed by the 
factor e 1 / 2 compared with the derivative with respect to y; 

(A) 2 -,. (20) 

Hence, the leading order equation of the order O(e ) contains only the y-derivative of the 0-th 
order solution (i.e., ot" Q y ■ ■), and 1st order equations are consisted of i-derivatives of 0-th 

order functions and y-derivatives of 1st order functions (i.e., ci(o), /%)),•'' > a (i)> ^(l)'' - ') as 
shown in Appendix ^] At each order, Einstein equations give ordinary differential equations for 
(am, • • • ) with respect to y. Then it is possible to solve the equations analytically. 

2.3 Solution at Zeroth Order 



The concrete form of the 0-th order equations are given in Appendix A-2 and the solutions at the 
0-th order are given as follows, 

<XQS)(v,t) = k#(y. *) + «(*)> (21) 

P(o)(y,t) = ^-^lnH(y,t), (22) 

7(0) (y,t) = lnd(t) + V2b$(t), (23) 

(o) (y,t) = -^\nH(y,t) + 4>(t), (24) 

where, 

g(y,t)^l- ^ + 8) d(t)a\y\, (25) 
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and 



'{Ir I) (20) 



Here d(t) can be identified with the scale factor on a positive tension brane and 4>(t) describes the 
behavior of the scalar field. d(t) describes the time evolution of inter-brane distance. There does 
not exist (3{t) because /3( ) corresponds to the lapse function, so we can always gauge away j3{t) by 
the redefinition of the time variable t. These zeroth-order solutions are quasi-static and their time 
evolutions will be determined by the first-order solutions. If we take a(t) = <p(t) = and d(t) = 1, 
we can retrieve the vacuum solution obtained in |52j-|54j. 

The above solutions are valid except for RS model with b = 0(A = —8/3). For RS model, the 
0-th order solutions are given by 

<*(f))(v,t) = -d(t)k\y\+a(t), (27) 

/3 (0 )(y,t) = -d(t)k\y\, (28) 

7{0) (y,t) = lnd(i), (29) 

<t>(o)(v,t) = 4>V), (3°) 

where k = \^2a/6 is the curvature scale of the bulk Anti- de Sitter spacetime. In the following, we 
assume 6^0 and RS model is treated separately in Appendix C. 
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3 Solutions for bulk metric and scalar field 

3.1 Solutions at First Order 
3.1.1 Equations 



The concrete form of the 1-th order equations are given in Appendix A-3. Substituting the solutions 
at 0-th order into the equations at 1st order, we get 



~(3«' (1) + ~ v^' (1) ) + I (V - I) (7(i) - v^ (1) ) 



= d 2 e 2V2b H H'^ 1 



3 J H 2 

2 , li2 _^°(J + 3d & + ^2 bd ^ 



a + 2d z + - 
6 



6 ir 



1 - 36 2 (J 2 ; 2 2 



+ 6 F ' (31) 



(!) " f - + v^' (1) ) + 1 (V - |) ^( 7(1) - v^ (1) ) 



d 2 e 2y/2 ~ b H 



St 2 + V2&6 - -4> 2 + -k 
6 a 



\/2<r / d 2 

1 - 36 2 a 2 - 2 2 



+ qV , (32) 



2a (i) + ^(i) " !T^ 8a (i) + 4/3 (D " 3 ^(D + 3 ^(i)) 



+ ( 6? -I)^(7(i)-^(D) 

2d + 2V2b&$ + 3d 2 + \/260 + f 26 2 + - J 2 + 2^d + ^ + 2V2b^ 

\ 2 1 d d d 



V2 
6 



Xf^y f 2d + 4\/26# + 4dd + 



d 



1 - 36 2 a 2 - 2 2 



+ 2^ > ( 33 ) 



11 



a' (1) \^4 + ~ d 



= «(1) + 



V2ad 



jj (7(1) - yfi-H 



(i) 



v/2 cr 

Th 



dy + &\ - q' (1) ) , (34) 



„ cr<2 



2\/2 

6(3a' (1) + - 7(i)) + -3-^(1) 



f 2^26 (V-f)^( 7(1) -V260 (1) ) 



Pf 2 



d 2 e 2 



+ -3-6(1- 36 --- 



• (35) 



These equations can be regarded as the ordinary differential equations for (am, /3m, 7(1), 0m) where 
the right-hand sides of the equations act as sources. Then the solutions are given by homogeneous 
solutions and particular solutions which are determined by d, (j), d and V(0(o))- 

3.1.2 Homogeneous Solutions 



The homogeneous solutions are given by 



a 



(i) 



b 2 - 2 

K a (t)H^ +F a (t) + P(y,t), 



2 2 

~3 



6^- 



/3 (1) = Kp(t)H^-+F p (t)+P(y,t), 
7(1) = K^H-^ + >/2M^(t) + 65 2 P(y, t), -3a/2^ P'(y, i), 



0(i) 



,,2 2 
_3 



K^HS^ + i^(t) + 3^6 P(j/,t). 



(36) 
(37) 
(38) 

(39) 



Here, F a , Fp, F^ and P are arbitrary functions and they do not affect the resultant effective 
equations. On the other hand, K a , Kp, Kj, carry the information of the bulk. These functions 
satisfy two constraints; 



3K a + K I3 + K~ / - 2^2b Ka 



0. 

2 

3A 



(40) 
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These constraints are insufficient to determine the behavior of the homogeneous solutions. The 
homogeneous solutions will be determined by the junctions conditions at the branes. 

We should note that in the case of RS model (b = 0) without scalar field = 0, these 
constraints are sufficient to obtain the concrete form of the homogeneous solutions. We obtain the 
equation for K as 

\K a + ~if 7 J + \K a + IkJ) = 0. (41) 

Thus we find 

K a + -K^ = l -Ce~ 4 \ (42) 

where C is the integration constant. It is well known that C is related to the mass of the Black 
hole in the bulk and this term will induce "dark radiation" on the brane. 

3.1.3 Particular Solutions 

From eqs. (|31[l - (j35jl . the particular solutions is given by 

a (1) = f a {t)H 2 -^ b2 \y, t) + g a (t)H (y, t)y + K^H' 1 '^ (y, t)y 2 + k a (y, t), (43) 

= f p {t)H 2 -^ b2 \y, t) + gp (t)H ^/a*) (y, t)y + hp(t)H~ l/ ^ (v, t)y 2 + kp(y, t), (44) 

7(i) = f,(t)H 2 -^ b2 \y, t) + g^t)H l -W#)(y, t)y + h^H' 1 ^ (y, t)y 2 + yfibk^y, t),(45) 

= /^t^"^^,^^^ (46) 

The functions (f a , g a , h a , fp, gp, hp, / 7 , g 7 , h 7 , fy, g^, h^) are determined by a(t), cf)(t) and d(t). 
The solutions are presented in Appendix B. On the other hand, the solution k(y, t) is determined 
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by bulk potential V as 

k' a (y,t) = l^^m H -2/^ f V dyH 2 ^V, (47) 
6 Jo 

k' p (y,t) = -Itfe^H-*' 3 * £ dyH 2 ' 3 * (V^ + (1 + 26 2 )u) 

_ ^ a - l de 2 ^HV, (48) 
k'^y,t) = -^d 2 e 2 ^H~ 2 /^ j\yH 2 ^(^ + V2bvY (49) 

3.2 Derivation of effective Priedmann equation 

Here we will demonstrate how we can derive the Friedmann equation using the above solutions. 
We use the solutions which we have derived in the previous subsection and the junction conditions 
given in Appendix lAl 



Projecting a',-^ onto y = 0, we get 
c/ (1) (0,t) = c$j = -yV<*(l - 6b 2 )f a {t) + g a (t) + k' a (0,t) - ^ad{3b 2 - 2)K a (t) + P'(0, t). (50) 
On the other hand, the junction conditions l)121j) gives 



(l) „, 2s , ,, s . . ,„ , x \/2 



°S = " V26^(t) + / 7 (t) - V2U(t)) + P'(0,i) - ^(p + U)de^. (51) 



So we notice 



2 



^-<rd(l - 6b 2 )f a (t) +g a (t) + ^-(Mt) - >/26/*(t)) + -(p + U)de^ 

= -— ad[-2(36 2 - 2)KM) + KJt) - yftbKAt)}. (52) 
6 

There appears the homogeneous solution in ea.l|52|). so we have to know the explicit form of it to 
obtain the effective Friedmann equation. In order to do so, we use the junction condition at y = r. 
Using it, we can write down the homogeneous solution in terms of the quantities on the negative 
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tension brane. Indeed, projecting a',^ onto y = r we find 



«(i)M) = «(i) 
"V2 



<rd(l-6& J )/ a (*) 



H 



l-(l/3ft 2 ) 



+ 



crd(l-36 2 )c/ Q (t) + 2/i a (t) 



H 



-1/362 



+ ^adh a {t)H~ l ^ b \ 2 + k' a (r,t) - ^ad(3b 2 - 2)K a (t)H- 2 '^ + P'(r,t), (53) 



and the junction condition ()124|) gives 



(i) 
(i) 



■^[(K 7 (i) - V2bK (f) (t))H- 2 / 3b2 + (/ 7 (t) - V2bU{t))H l -^ b2 ) 
+ ( 57 (t) - V26^W)^- 1/3b V + (h y (t) - V^bh^H- 1 ^ 1 



+ P'(r,t) + '-^-(p + U)de V2b t (54) 

D 



Combining them each other, we obtain 



^ad(l - bb 2 )f a (t) + g a (t) + ^ad(/ 7 (t) - v^/^i)) 



l-(l/3fe 2 ) 



+ 



\/2 a/2 
^-ad(l - 3b 2 )g a (t) + 2h a (t) + ^-ad{ gi {t) - yfibg^t)) 
3 b 



#-1/36% 
#-l-(l/3fe 2 ) r 2 



\/2 \/2 
+ ^—adh a {t) + ^—ad{hJt) - V2bhJt)) 
3 b 

+ k' a (r,t)-^(p + U)de^ 
b 

= -^ad[-2(3b 2 - 2)K a (t) + K 7 (t) - \/2W^)]i^ 2/3fe2 . (55) 

Explicit form of all homogeneous solutions are given in Appendix^! Now let us use the homogeneous 
solution (|55|) . Substituting (|55j) into (|52j). we finally obtain the effective Friedmann equation. The 
concrete form of the effective Friedmann equation (|58|) is given in the next section. Similarly, we 
can derive other effective equations. All of them are also given in the next section. 



Here we note that one can not determine the homogeneous solution in one brane model because 
one knows the only one junction condition at the brane where one lives. 



15 



4 4D Effective equations 



4.1 Priedmann equation, equations of motion for radion and scalar field 



Now we can explicitly write down four-dimensional effective equations. First, the evolution equa- 
tions for scale factor are given by 

26? + U 2 = ^V((l - 3w)p + 4U + 3vW)e-^ - -W 0) ), 



" f '-" ' 6" 36 " ~ vv ~ ~~ /r ' " ' 6 

d + 26? + -<^ 2 = -— K 2 a((l - 3w)p + 4U + Sy/Mj^e-^ H~*^ 



(56) 



36 



« 2 + ^-^= ^ (A t 4 ^ V 



6" 24(1 -HC) 
V2 



3A + 8 

p + U + (p + U)H3K+8 



(3A + 8) 2 



(57) 



-V2b<}> 



— (A + 4W ^- /" dyH^V 

24 V ' i- H ( Jo 



+ 



3ff^(A + 4) 
1-tfC 



1 



3A + 8 



(3A + 8) 2 



d 2 H~ 2 



The equation for (p is given by 

bn 2 o 



4>+3&0+V2b0 2 



6(1 - #<) 
V2 



+ 



(l-3u;)p+4[/-v / 2r 1 f/ / +((l-3tZ')p+4C/-v / 26~ 1 L> / )^ 3 ^ 

1 /*^* 8 

A + 4)ad- — J dyH^K (y' + y/2bV) 



(58) 



-V2b4> 



1 v^b 



3A+16 V 



+ 



l-Ht 



1-H< 2 
6V2b 



V^) 1 -#W8 



(0 (O) A 

^(o) )y | 



3A + 8 



3A + 8 



(3A + 8) 5 



(59) 



Finally the equation for d(t) is given by 

d + 3d& + \/2&# + O ? ^ + 4 v d 2 g~ 1 
z(oA + 8) 

- ^ 3A + 8 1^ K V [(l _ 3w)p + 4£7 + 3^26*7' + ((1 - 3w)p + 4U + 3V2bU')H~^ 



72 



8 + 3A 
12~ 



HV{(f>( 0) ) ( 1 - iJs+k 



^(o))^ 



, (60) 
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where we defined 



H = H(r,t) = l-d(t), 

3(A + 4) 
^ 3A + 8 ' 

To make the spacetime non-singular, we set the negative tension brane at r 
assume d < 1. We can also get the energy conservation law as follows, 



p + 3(1 + w)p a 



p + 3(1 + w)pa = 0, 

^ M- 1 ) = o, 



(61) 
(62) 

l/(V2b 2 a) and 

(63) 
(64) 



3A + 8 

which can be derived from above equations. Because we have five equations for five unknown 
functions (a, (j), d, p, p) we can consistently solve these equations. 



In the RS model, the effective equations are given by 



& + 2d 2 + x ; ;J 



d + 2d 2 + 



6 r 
1 v 



a 



2 1 22 



+ 3d</> 
(i + 3dc? — d 2 
Here we take r = 1/k. 



6 -((l-3w)p + 4U) r 
-^((1 - 3*)p + 4L>) + d + 3<id - (i 2 - e" 2 ^, 



k 2 /c e 



2d 



3 e 2d -l 



p + f/+(p+?7)e- 4d +— V(l + e 



1 



-2dN 



3 2rf 



1 



-K 2 k- 

( 

K 2 fc r 



(U' + U'e- M ) - -(l + e- 2d )V 



,2d_ 1 



2d 



((1 - 3w)p + AU + ((1 - 3w)p + 4C/)e 



-2d 



(65) 
(66) 

(2M-d 2 ), (67) 

(68) 
(69) 



1 

^(l- e - M ). 



4.2 Effective Theory 



In the previous section, we derived effective four-dimensional equations. We try to interpret them 
in terms of the scalar-tensor gravity. In this section the effective theory for observers confined on a 
positive tension brane is considered. In this section, we do not consider the potentials U, U and V. 
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The action for the scalar-tensor theory is given by 



<P 



(70) 



The Friedmann equation and the equation for cp are given by 



• 2 ,if> (j(<p) ( if 

a + a— 



(p + <p 3a + 



6 



2w + 3 



3 ip ' 



2uj + 3 



;i - 3w)/9. 



(71) 



4.2.1 Casel: d = const. 

At first, let us consider the case where the radion is somehow stabilized, that is, d = 0. In such a 
case, we find that eq. becomes 



a 2 + v 7 ^ 



n c, 



1 

6 ( 



+ 3A<£ + \/260 2 



V2(A + 4) 
24(1 - fft) 

6k 2 ct 
6(1 - #f) L 



(1 - 3w)p + (1 - 3w)pH3S+s 



If we take 



(72) 
(73) 

(74) 



3A + 8 2b 2 ' 

the equations become the same as (|71|) . Then the effective theory can be identified with Brans-Dicke 
theory with Brans-Dicke parameter oj given in (|74|) . We also find the four-dimensional gravitational 
constant becomes 



>/2(A + 4) 



K 2 a, K 4 = K 2 i?3A+8. 



(75) 



* 8(1 - H<) 

Here H is constant because d is constant now. It is interesting that the matter on the negative 
tension brane couples to the gravity with a different gravitational constant from the matter on a 
positive tension brane. 

From the observational constraint, we need co > 3000 at least at the late stage of cosmology 
|55| . So the dilaton coupling is strongly constraint as b 2 < 1.6 x 10 -4 . Thus in the HW theory with 
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6 = 1, the stabilization mechanism for eft should be introduced in order to obtain an acceptable late 
time cosmology. 

In RS model, the theory becomes Einstein theory and the four-dimensional gravitational con- 
stant is given by a well known formula; 

2 2 e 2d 2 2 4d 

4.2.2 Case2: = const. 

Next, we consider the case where the bulk scalar field is stabilized, i.e., if> = 0. This case has been 
studied by many authors in the context of the Randall- Sundrum model (20], [SJ, In such a 

case, ea.(|58|) becomes 



, 2 ff<(A + 4) 
a — 6 = 

i -Hi 



1 -MH- l + - - 1 S J 2 H- 2 



3A + 8 (3A + 8) 2 



V2(A + 4) c _^ b$ 
24(1 - H<) 



p + pH s a+s ; (77) 



J + 3<M + x/26# + 3Z ^ + 4 d^-i = V2(3A + 8) ^ 2(7 r (1 _ ^ + (1 _ 3t2i) -£-§tgl e -V&* (78) 
2(3A + 8) 72 L 

From these equations, if we take 

— ^ « = STl(l^)=WTij(i^)- < 79 » 

eas. (|77|) and (|75|) can be written as (|T1|) . Now w is a dynamical variable, so we conclude that 
scalar-tensor gravity is realized. The effective four-dimensional gravitational constant becomes 

w 2_ V2(A + 4) ^ R 2^ K 2 il _ (p)J( £^ = K 2 {1 _ lp)1 J ri _ m 

o 

From the observational constraint uj > 3000, we find 

( 1 \ VC 



1500(A + 4) + 1 

The physical distance .D between the positive and negative tension brane is given by. 

i 



D = dyd e v ^ = -p— e VZb<f >. (82) 
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Thus the constraint means the negative tension brane should be sufficiently far away from positive 
tension brane. For b 7^ 0, the position of the singularity corresponds to d = 1. Thus the constraint 
requires that the negative tension brane is close to the singularity. 

In RS model with A = —8/3, the BD field and BD parameter are given by 

v= i-.-», -l(^) <«> 

and the four dimensional gravitational constant is defined by 

K\ = n 2 k, k\ = k\{1 - iff. (84) 

These results completely agree with the results obtained in |5Uj.|51j.|56j. The observational con- 
straint implies 

D = d/k> 4/k. (85) 

4.2.3 General Case 

In general cases, the effective theory becomes bi-scalar tensor theory. For small b satisfying the 
constraint b 2 < 1.6 x 10~ 4 , we can obtain a phenomenologically acceptable cosmology without 
introducing a stabilization mechanism if the negative tension brane is sufficiently away from the 
positive tension brane. However, the dynamics of the moduli fields strongly depends on the matter 
contents on both branes and the initial conditions for the moduli fields. Thus in many cases, it 
seems to be inevitable to introduce some stabilization mechanism for moduli fields. The stabi- 
lization problem of moduli fields are very general and serious problem which appears in all higher 
dimensional theories. Detailed analysis of the general dynamics of the moduli and investigation of 
the stabilization problem is beyond the scope of this paper. We will come back to this issue in 
future publications 57 . 
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5 Boundary inflation 



In this section, we consider a simple model for inflating branes driven by potentials U, U and V. 
There are essentially two types of the inflation [31 . First, the potential energy of the brane potential 
U can drive the inflation on our brane. This can be called "matter field inflation". Second, the 
potential energy in the bulk can also drive the inflation on our brane. This can be called "modular 
inflation" . In general, one could have a mixture of both types of the inflation. In the following, we 
will provide a simple model for "matter field inflation" and "modular inflation" respectively. For 
simplicity, we do not consider the matter on both branes, that is p = p = through this section. 

5.1 Matter field inflation 

First, we consider a simple model for inflation driven by brane potential U. Then we assume that 
the bulk potential vanishes V = 0. Furthermore, in order to avoid the negative tension brane hits 
the positive tension brane or singularity, we assume the moduli fields d(t) and cj)(t) are stabilized 
by tuning potentials U and U. From the equations of motion for d(t) and (j)(t), we notice that this 
can be achieved by choosing 

U = e 2 ^U , U = -UhI£tI, (86) 

where Uq is independent of <f>. Then it is possible to stabilize the moduli fields d(t) and 4>(t); 

d(t) = d* = const. 4>{t) =(/>* = const. (87) 

The effective equations become a simple equation 

& = ^±Q e&+-K?*U . (88) 

If one considers that the potential energy Uq is provided by some matter field on the brane, the 
conventional four-dimensional inflationary scenario can be realized. We comment on the stability 
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of the moduli fields d(t) and 4>(t). Let us consider small fluctuations of the moduli fields; 

d(t) = d* + Sd(t), <f>(t) = </>* + 5<f>(t). (89) 
Linearizing the equations of motion for 4>{t) and d(t), we get the equations for fluctuations; 

Sd + 3&6d + (3A + 4)6> 2 Sd = 0, (90) 

• • 6\/2 H** ■■ ■ ■ 

5<f> + 3a5<f> = ~ 3A + 8l _ jy C 0* d + 3&Sd )- ( 91 ) 

Then we found that for 3A + 4 > 0, the fluctuations acquire a positive mass squared, thus the 
solutions are stable. On the other hand for 3A + 4 < 0, the system of inflating two branes is 
unstable. This agrees with the result that in RS model with 3A + 4 = —4, the inflating two branes 
are unstable and the radion fluctuation 6d(t) has a mass squared —4a 2 . It would be interesting to 
quantify the effect of these moduli fluctuations on the inflationary scenario. 

In RS model with no bulk scalar field, it is possible to construct a viable model without the 
stabilization of the moduli field d(t). Let us consider the matter field inflation driven by potentials 
on a positive tension brane U. The effective equations are given by 

d + 3ad-d 2 = U. (93) 



The solution for U = const, can be easily found as 

a 2 = d 2 = ^U. (94) 

The fate of branes depend on initial conditions. If the initial condition for d is properly chosen, the 
inflation on positive tension brane pushes the negative tension brane away. Thus the observational 
constraint on d becomes easy to be satisfied. It is also interesting to consider a model where branes 
will collide, although there are many open questions in cosmology based on colliding branes. 
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5.2 Modular inflation 



Next let us consider an inflation driven by bulk potential V. In this case, the moduli <fi acts as an 
inflaton. Thus we should not stabilize <j>. The stabilization mechanism should be switched on after 
inflation. In general, the other moduli d(t) also participates in the dynamics of inflation. In order 
to simplify the model, we take a bulk potential as 

V = -Ae-z^o 2 . (95) 

Furthermore, as is done in the previous subsection, we stabilize d(t) by tuning the brane potential 

U on the negative tension brane as 

~ = 2_V2aA - m ,_ _ v - 2b ^ 
fir A 

Then we can have d(t) = d* = const, and the effective equations become 

& + V2b&i-U* = -^-^aAe~ 2 ^ b i (97) 
6 6A 

} + Mi + V2bj> 2 = -^aAe- 2 ^ b t (98) 

This equations yield a power-law solutions for scale factor e a and scalar field e^. An interesting 
property of these solutions is that the power does not depend on the amplitude of the potential 
A and it is determined only by the dilaton coupling b. We find that a power-law inflation can be 
realized for —2 > A. 

Finally, let us consider a RS model. In this model, we can simply take d — > oo to neglect the 
contribution of the moduli d(t) rather than stabilizing d(t). Then the effective equations become 

6^ 12 

+ 3&<f> + -V = 0. (99) 

If we define 04 = 4>/ Kk 1 / 2 and V4 = V r /4« 2 fe, these are completely the same as the equations obtained 
in conventional four-dimensional theory. Thus the bulk inflaton can mimic the four-dimensional 
inflaton dynamics 11 . 
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6 Conclusions and Future Works 



In this paper, we considered two branes model with bulk scalar field, which is a generalization of 
both Randall-Sundrum model and Hofava-Witten model. The bulk potential of the scalar field has 
a exponential functions Vb u ik = exp(— 2y/2b<p). We considered matter on both branes and arbitrary 
potentials in the bulk and on the brane. These matter and potentials induce the cosmological 
expansion of the brane as well as the time evolution of the bulk scalar field and radion. Because it 
is expected that the tension of the brane is much larger than the energy density of matter at late 
time cosmology, we have used the low-energy approximation to derive the equations which describe 
the cosmological evolution on the brane. 

The main results of this paper are the four-dimensional effective equations (|56|) - (|6U|) which 
govern the low energy dynamics of the two branes system. The important point is that once the four 
dimensional dynamics is determined by these effective equations, we can construct a correspondent 
full five-dimensional geometry. This correspondence would be useful to investigate the effect of the 
bulk geometry on four dimensional brane dynamics. 

We also investigated observational constraints imposed on the model. We assume we are living 
on positive tension brane. First, we studied the case where the radion d(t) is stabilized. In such a 
case, we can regard this model as Brans-Dicke theory with a Brans-Dicke parameter 1/26 2 . From 
the observational constraint, we concluded b 2 must be smaller than 1.6 x 10~ 4 . Second, we consider 
the case where 4> = const. We found that this model can be identified with the scalar-tensor theory 
where the Brans-Dicke parameter is given by (3/2(3b 2 + l))(ip/l-ip), where ip is the BD field which 
is defined by d(t). For RS model with 6 = 0, this agrees with the result obtained in |5U] . |5T] . [55] . 
The observational constraint implies the negative tension brane should be sufficiently away from 
the positive tension brane. In general cases where there in no stabilization mechanism, the effective 
theory becomes bi-scalar tensor theory. 
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Using the effective equations, we constructed severai simple models for inflation. The inflation 
can be driven by potentials on the brane and bulk. We can construct various models for bound- 
ary inflation. It is very interesting to know we can distinguish these models from conventional 
four-dimensional model. The detailed analysis of these inflation models including the primordial 
fluctuations will be given in the near future. 

Note added 

After almost completing this work, we found jHS] . In jHH] , the equations for the radion and the bulk 
scalar are derived using a moduli approximation and qualitatively the same result is obtained. 
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Appendix 

A Five-dimensional equations 

A.l Einstein equations and scalar field equation 

Five-dimensional Einstein equation and scalar field equations are given by 



(y,y) ■ 



a' 2 + a'0 - (a + 2d 2 - a$)e 2 ^ 



6 



1 



0'2 + I^ e 2( 7 -/3) 



I. 



1 



(■ 



6 2 _ A e 2( 7 -V2^) a 2 _ V^yt] . ( 100 ) 
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(M): 



a" + 2a' 2 - 7V - (d 2 + -yd)e 



2(7-/3) 



6 6 6 



+ \e> {[-y/2e~^a - K 2 {p + U)]S(y) + [V^e'^a - K 2 (~p + U)]S(y - r)}. (101) 



2a" - 2j'a' + 2a'/?' + 3a' 2 + 0' + /?' 2 - P'j' 

- (2a + 2 7 d - 2a$ + 3d 2 + 7 + 7 2 - /?7)e 2{7_/3) 



.V 1^2 2( 7 -/3)_ 1 

V V 2 



6 2_f ) e 2(7-^) CJ 2 _ y((/))e 2 7 



+ e 7 { [-v^e-^<T + « 2 (p - [/)]<%) + [V^e'^a + n 2 {p- U)]S{y - r)}. (102) 



(y,t): 



a' - ja' - a((3' - a') = -\m>' ■ 



(103) 



Equation of motion for 



b" + 3a' 4>' + p'cj)' - 7>' - {4> + 3d<^> + 70)e 2 ( 7 -^ 

= _ ^26 f 6 2 - -) e^-^a 2 + i^e 27 
\ 3 / 2 a<p 



+ e 7 



-2be~^a + k 2 ^-) 5(y) + ( 2be~^o + k 2 ^ | <5(y - r) 



(104) 



We expand these equations in terms of e. 



A. 2 Oth-oder 



The field equations of the 0-th order are given as follows, 



«(0) + a (0)P{o) - -M) - fi 



5 Z _ Z I cr 2 2(7(0,-^0(0)) 

6 V 3.' 



(105) 
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«8» + Ho) ~ ioAo) = - I (&* - I) ^co)-^(o)), (106) 

2a ('o) - 2 7( ) a (o) + 2a (o)/ 3 (o) + 3a (o) + tf[o) + ^(o) - ^(0)7(0) 

= -^f 0) - \ (b 2 - I) a 2 e 2 ^(o)-^(o)), (107) 



»(o) - 7(o)«(o) ~ "(o)(/?( ) - a (0)) = -o ^(o)^(o) > ( 108 ) 



4>'[ 0) + 3a' (0 / (0) + (3{ 0)( f>' m - 7 [ O )0(o) = -^26 (&* " |) a 2 e 2 ^(o)"^(o)). (109) 
Next, we can get the junction conditions at y = as follows, 

a §)W = -^^(°)-^(o)| y=0 , (HO) 

^(O)^) = -^^(0) -^(Ol^o, (111) 

$>(<) = -6a e^-^O) 1^=0. (112) 

Similarly, we can get the junction conditions at y = r as follows, 

= -^^(o)-^(«)|, =r , (113) 

= -^ae^)^^)|, =r , (114) 



'(ojv; - 6 



#JJ(t) = -ba e^)-^(o)\ y=r . (115) 
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A. 3 1st order 



In this subsection, we consider the equations of the order e 1 , which determine the behavior of 
a(t),(f>(t) and d(t). They are given as, 



2a (0) a (i) + a {o)P(i) + a (i)^(o) 



ci(o) + 2d 



(o) 



d (o)Ao)J e 2 ^)-P m ) 

~l(b 2 - |) ^ ( 7(1) - V2b4> {1) ) e 2 «°)-^(°)) + i e 2 W, 

(116) 



«(i) + 4 «(o)«(i) - 7(o) «(i) - 7(i) «(o) " («(o) + 7(0) "(0)) e 2(7 <°) 

= -^(o/(D " ^o)^ (7(0) -* 0)) " I ( b2 ~ I) (7(D " V^ (1) ) e 2 ^(o)-^ ( o)) + l e a7 ( o,v; 

(117) 



2 «(i) " 2 7( 0) «(i) " 2 ii) a '(0) + Ho^i) + 2 «(i/(o) + 6 «(0)«(i) + ^i) + 2 P{oAi) ~ ^(0)7(1) - /3(i)7( 0) 
+ (-2a (0) - 27 (0) d (0) + 2d (0) /3(o) - 3d 2 0) - 7(0) - 7( 2 0) + /? (0) 7(o)) e 2 ^)"^)) 

= -0(o)0(i) + ^?o)« 2(7(0) ^ (0)) " {f " |) - 2 (7(1) " >/2ty(i>) e 2 ^o)-^(«o) + I e 2 W, (118) 



«(i) - (7(0) «(i) +7(1)0(0)) 



«(o)(/3(i) - a' (1) ) + d(i)(/3( ) - a(o)) = - 3 (0(o)0(i) + 0(i)0(o)) > ( 119 ) 



0(i) + 3a (o)0(i) + 3c *(i)0(o) + ^(0)0(1) + ^(i)0(o) - 7(0)0(1) - 7(i)0(o) 

- (0(0) + 3d( O )0(o) - /3(o)0(o) + 7(0)0(0)) e 2(7(0)_ ^°) ) 



-2y/2b (b 2 - - ) a~ [-y (1) 



V2b4> {1) ) e 2(7(0)-^ (0) ) + 1 ^ (12Q) 
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Similarly, we can derive junction conditions at y = and y = r as follows, 

a (X) W = ""T ( 7(1) " ffe7(O) ^ fe0(O) " y (P + ^) e7(0) ly=0> ( 121 ) 

= -^(7(i)-v / 260( 1 ))ae^)-^(o) + ^(3p + 2p-C/)e^) | y=0 , (122) 

= -K^)-^a)) ffe7w "^* w + T^ 67(0) U' (123) 

ag(t) = -^( 7 (i)-v / 260 (1) ) ( 7e^)-^(o) + ^(p + C/) e 7(o) ^ (124) 

= -^(7(i)-v / 260 (1) )^(o)-^(o)-^(3p + 2p-[/) e ^) 1^, (125) 

$>(t) = -6(7(i)-^(i))^>-^)-y^e^) \ y=r . (126) 
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B Solutions for 1st Order Equations 



B.l Homogeneous solutions 



Here we will write down the homogeneous solutions in terms of p, U. 



.2^-1 



K n = - 



3A + 8 
2V2a- 2 



<P + U) 



,2V2b<j> 



A 

6^2 



,V2b4> 



, _ 3(A+4) 
- A 1 ff 3A+8 



3A + 8 



daH^+s + 



(3A + 8)' 



-„ _ 3A+4 
-d H 3A+8 



K R = 



V2 



-2^-1 



*((3.o + 2)/!-C;)e^ 1 * 



3A + 8 



A 

6^2 / d 



3A + 13 ~ 2 \ .-. 3(A+4) 

^ ) # 3A+8 



3A + 

A 
A 



- 4 6\/2 ; 9 - 3A+4 

.^if 3A+8 



- + 2y/2b-4> - , 4 

d dW (3A + 8 



(127) 



r -—^H 3A+s / dyH^+s 



V2bV' + (1 + 2fe 2 )y 



(J 



A 



e 2 ^HV, (128) 



(129) 



A 

6\/2 



f + 3d0 + v^ 2 - 3^26 (& + 2d 2 + \i 12 
(A + 4) V V 6 



366 



3A+4 



H 7 d0^3A+8 _ -77d 2 i? 3A+8 

3A + 8 Y (3A + 8) 2 

_^/2^d e 2^ # -3^3 r d yH*&*(y + y/2W). 
A Jo 



(130) 
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B.2 Particular solutions 



The particular solutions for 1st order equations are given in eas. (|43|) - (|46|) . The concrete form of 
them becomes as follows; 



fa 



n 

9a 



K 

hp 
/?.,. 



4a~ 



(A + 4)(A + 2) 

4a~ 2 
(A + 4)(A + 2)' 



4a~ 2 
(A + 4)(A + 2) 



^2V2b4> 



a 2 + \f2bac 



1 22 3(A + 4) / d x Id 



2a' + a 2 + y/2b$ + 



- a J 

6 3A + 4 \^d 2d 2 

3A + 13 79 



6 



3(A + 4) 
3A + 4 



/-, d i d Id 2 



^2V2b4> 



+ 



3&<j> + V2b<t> 2 - 3V2b (a + 2a 2 + ^<p 2 ^j 



3(A + 4) (d j 3^/2, d 2 " 



U = V2bL 



3\/2a 1 2 V2bj> 
' 3A + 4 



2da + — 



D 2V2bj> 



3A + 4 

18a- 1 
~3A + 4 f 

V2bgs, 



_1 J2 „2V2b^ 

4 a e 
4 a e 



e 2V260 [2d + 4V2bd4> 



4 + f4 



4 



(131) 



, (132) 



, (133) 
(134) 
(135) 

(136) 

(137) 
(138) 
(139) 
(140) 

(141) 
(142) 



For A = — 2 and/or A = —4/3, the particular solutions should be treated separately. But even 
in such cases, we get the same effective 4D equations as we have shown in Sec. 4. 
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C Homogeneous and Particular Solutions in RS model 

C.l Homogeneous solutions 



P(i){y, t) 
7(i) (y, t) 
<t>(i){y, t) 



K a (t)e 4kd y + F a (t)+P(y, t), 
Kp{ty kd y + F p (t) + P{y, t), 



K 7 (t)e 



Akdy _ n/ 



d 



P'(y, t), 



= K^e 4 ^ + F <p (t), 



(143) 
(144) 
(145) 
(146) 



where F a ,Fp,F^ and P are arbitrary functions which do not affect the resultant 4D equations. 
K a , Kp, K 1 and K$ yield the the constraints; 



3K a + Kp + K 7 = 0, 



(147) 



K a + (K a - Kp)& + ^K^ + = 0. 



Here we define k = (\/2/6)cr as in Sec. 2. 3. 



C.2 Particular solutions 



(148) 



The particular solutions are given by 



0(1) 
7(1) 



4fc 2 



• 2 1; d- Id 2 
a <p + -a + 



6 d 2d? 2k 



+ — ( da + 



Id 2 



1 ■ 



y - V 



2 d 2 / a I 
id 2 \ i 



2fcdj/ 



I ■■ X9 5 a d 1 </ 2 \ 1 / •• . . «, , 



(149) 



e 2My , (150) 

(151) 
(152) 
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